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ABSTRACT 

The  narrowband  ambiguity  function  has  been  used  to 
study  the  space-time  resolution  of  radar  and  sonar  signals. 
Its  derivation,  properties,  and  generalizations  are  reviewed 
in  a  unifying-  notation.  A  "most  general  ambiguity  function” 
is  defined,  from  which  both  narrowband  and  wideband  am¬ 
biguity  functions  can  be  derived. 
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A  REVIEW  OF  NARROWBAND  AMBIGUITY  FUNCTIONS 


INTRODUCTION 

The  narrowband  ambiguity  function  was  introduced  by  Woodward  (1)  in  1953.  Its 
properties  have  been  extensively  explored  and  used  to  study  range  and  velocity  resolu¬ 
tion  of  radar  and  sonar  targets  by  Woodward  (1)  and  others  (2-41).  The  wideband  am¬ 
biguity  function  is  of  more  recent  origin  and  has  not  yet  been  as  fully  developed.  Refer¬ 
ences  to  both  functions  are  scattered  through  the  literature,  in  v, klo.ty  varying  notations. 
This  report  reviews  the  derivations  and  properties  of  the  narrowband  [unction  and  some 
of  its  generalizations,  using  a  consistent  notation.  A  companion  report  ("A  Review  of 
Wideband  Ambiguity  Functions,")  is  concerned  with  the  wideband  function  and  the  rela¬ 
tionship  between  the  two  functions. 

These  reports  originated  as  a  set  of  notes  for  a  Branch  Seminar  and  have  been  re¬ 
vised  for  more  widespread  distribution.  The  intent  was  to  elaborate  on  and  explain 
Woodward's  (1)  Chapter  7  and  then  to  discuss  some  of  the  subsequent  work.  To  facilitate 
reference  to  Woodward,  the  initial  sections  follow  his  work  closely.  Notational  changes 
were  necessary  for  consistency  in  the  sequel.  Details  of  some  of  the  developments  have 
been  relegated  to  the  appendices. 


RESOLUTION  AND  AMBIGUITY 
Range  Ambiguity 

Let  a  signal  transmitted  at  time  f  be  represented  by  the  real  part  of  the  analyticf 
signal  f(  t).  We  consider  first  the  problem  of  range  resolution  of  point  targets,  where 
range  is  determined  by  the  known  velocity  of  propagation  and  the  measured  delay  in  the 
signal  echo.  We  assume  no  attenuation.  To  achieve  maximum  resolution,  we  would  like 
the  echo  f  (  r  -  t)  to  differ  as  much  as  possible  from  f(  t).  Using  a  mean-square  cri¬ 
terion,  we  would  liket 


i  |f  (t  )  -  f((-T)|J  dt  (1) 

to  be  as  large  as  possible,  except,  of  course,  near  r  o.  That  is,  we  wish  to  maximize 

/[f  (t)  f(  f  -  T)]  [f*(t)  -  **(«-  T)]  dt 
=  f  !f  (  t  )j  2  dt  *  j',f  (  t  -  T  )  |3  dt  -  f<V(t)'V*(t  T)  dt  ~  /  f  (  f  -  T)  f*(()  dt 
IE  2Ki  [I  f  (  t)  f*(t  -  t)  o'(] 

where  i|f(()|3cit  ljf(<-T)|3  dt  t,  the  total  energy  of  the  analytic  signal  or  twice 

the  total  energy  of  the  real  signal.  Equivalently,  we  wish  to  minimize 


+A  complex  signal  is  said  to  be  analytic  if  its  imaginary  part  is  the  Hilbert  transform  of  Its  real 
part:  in  [f  (  t )]  -  «  <sf  [f  i  « )] )  . 

I  Here  and  In  the  sequel  integration  Is  over  the  entire  space  of  the  appropriate  dimension,  ±»  for 
all  single  integrals,  unless  otherwise  stated. 
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*  •  Iff  (t)’l’V--r)  rfO  •  (2) 

except  near  *  =  o.  Let 

*(*)  =  “(O  exP  (iut)  .  (3) 

Then 

»*(t-T)  =  U*(t  -  T)  exp  -  T)]  . 

and  Eq,  (2)  becomes 

B-t  [exp  (i«T)  /u  ( t )  u*(  t  -  T)  dt]  (4) 

which  oscillates  with  r. 

The  requirement  is  that  1* (r)|  be  as  small  as  possible,,  where 

*(r)  *  f  u(t)  u*{t  -  r)  dt  (5) 

is  the  complex  autocorrelation  function,  and  the  symbol  -  means  "is  equal  by  definition 
to."  Sincet 

^tr[u(<~T)l  -Ju(t-r)  exp  (-2  nift)  dt  -  exp  (-2  -ni  fr)  U(()  , 

Parseval’s  theorem  gives  us  from  Eq.  (5) 

*CO  «  /|®(/)|a  exp(2>ri/T)  dt  ;  t|^(Jr)|a3  •  (6) 

Furthermore,  Jf  (0 )  =  /iu(t)|s  dt  -  / 1 1/ v /") j 3  df  -  E. 

As  a  measure  of  total  signal  ambiguity,  Woodward  (1)  defines  the  time-resolution 
constant 

T*  £/|*0‘)l*  dr  =  ^-/ROr  df  .  (7) 

where  the  last  equality  again  uses  Parseval's  theorem. 


Ambiguity  in  Range  and  Velocity 

If  the  effect  of  moving  targets  is  assumed  to  be  adequately  approximated  by  a  simple 
shift  in  frequency,  we  can  define  in  an  analogous  fashion  in  terms  of  the  frequency  shift  a 
"frequency  autocorrelation  function" 

K(«J>)  $  df  ® 

=  /|U(0|2e*P<-2ni®f)  dt  =  yt*[)u(t)|3]  .  (9) 


tin  general,  5\,r[  f  (*)]  =  t(y)  symbolizes  the  Fourier  transform  t  ( y)  =  f  f  <*>  r«p  <- )  Jx. 
The  inverse  transform  will  be  symbolized  by  tf  [  s  ( v  >]  *  f  <*)  =  I  <  <y>  *»p  (  2»i»r)  «/y. 
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Since  If  [l/(  /-♦)]  -  exp  u(  t ). 

Without  loss  of  generality  we  can  let  £  -  l,  and  by  analogy  with  Eq.  (7)  define  the 
frequency  resolution  constant 

F  *  /|«<<ti|J  d*  -  / !  u  f  r )  1 4  dt  ,  (10) 

where  the  second  form  comes  from  Eq.  (9),  using  Parseval's  theorem. 

If  the  targets  are  at  different  ranges  and  are  moving  with  different  radial  velocities, 
we  need  a  combined  time  and  frequency  correlation  function  (two-dimensional  correla¬ 
tion  function) .  For  the  time  being  we  consider  only  the  narrowband  approximation: 

The  Doppler  effect  is  approximated  by  a  frequency  shift  <t>,  constant  across  the  signal 
bandwidth. 

The  echo  of  t)  is  thus  given  by  *(  t  -  ts  exp  [27M<t(  (  -  r)].  Again,  we  let  »( t)  = 
u(t)  cxp(iut)  z  U(t)  exp(2*i/>  >.  The  function  to  be  minimized  is  now,  by  analogy  with 
Eq.  (2), 

gi  [exp  (  2wj<t>r )  /  ¥(<)  **<  t  -  t)  exp (- 2ni«l )  dt] 

=  J?  i  {exp  [2  ni{f  ♦  <t>)  t  ]/«•(»)  u*(t  -  r)  exp  (~2rri<t>!  )  dt]  . 

We  require  the  modulus  of  the  combined  time  and  frequency  correlation  function 

X  (r  ,<t>)  =  /  u  ( t )  u*(  t  -  t  )  exp  (-2 rriOt)  dt  (11) 

to  be  as  small  as  possible,  except  near  |x(0,0)(  £  1. 

The  ambiguity  function  (Woodward  ambiguity  function,  narrowband  ambiguity  func¬ 
tion,  n-b  autoambiguity  function)  is  defined  as 

-*(’-.*)  =  U(T-®)|2  (12) 

Other  definitions  of  the  generalized  autocorrelation  function  (GACF)  which  lead  to 
the  same  ambiguity  function  appear  in  the  literature.  For  example,  let 

A  ,  j, 

X  -  J  U(l  T  2)  U*(  t  -  T  2)  exp  (-2w«t>f )  dt 

-  |o(1)u*((-t)  rxp  [-  2i» /<!>('/  t '2))  dij 

-  exp  (tti4t)  V(t,®)  , 

so  that 

lx  !x(t,4>)!2  . 

The  present  definition,  Eq.  (11),  has  some  useful  transformation  properties  which 
will  be  discussed  later. 

The  name  of  the  ambiguity  function  stems  from  the  fact  that  it  does  not  uniquely  de¬ 
termine  a  waveform.  For  example,  let  the  GACF  of  Eq.  (11)  be  tabled  \„<r  >  to  dis¬ 
tinguish  it  from  v,  i  t  .  mm  mi  i  )t-4i>  2-H>»)  dt.  where  \  x  /  >  m,  n .  Ambiguity  is 

apparent  if  we  consider  the  frequency- shifted  and  Imu --delayed  waveform 
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v(t)  -  u  ( t  -  (  )  exp  (  2-niOt )  , 

where  <  *  0  and  «  i  o,  Then 

=  \Xu(r,*)\2  . 

SOME  PROPERTIES  OF  THE  GENERALIZED  AUTOCORRELATION 
FUNCTION  AND  THE  WOODWARD  AMBIGUITY  FUNCTION 

X(T,0)  =  /u(  t)  u*(f  *  r)  dt  =  K(r)  (18) 

X  (0,4>)  *  /|u(t)|Jexp(-2Triit>t)  dt  -  *  (<t>)  (14) 

X(0,0)  =  J|u(t)|2  dt  -  E  s  1  .  (15) 

Theorem  1 

The  GACF  as  an  integral  In  the  frequency  domain  Is  given  as 

X(t,<6)  =  /{/*(/)  V(  ft-  <t>)  exp  (  2ni  fr)  df  .  (16) 

The  proof  of  theorems  1  through  4  will  be  found  In  Appendix  A. 

Theorem  2 

The  combined  time-frequency  resolution  constant  is  given  as 

/Jk(T.*)|*  drd*  =  1  .  (17) 

“The  effective  ‘area  of  ambiguity’  in  the  time-frequency  domain  is  independent  of  the 
transmitted  waveform  and  is  equal  to  unity.”  This  theorem  is  called  the  "Radar  Uncer¬ 
tainty  Principle"  by  Siebert  (3).  It  is  one  of  the  most  important  properties  of  the  narrow- 
band  ambiguity  function  and  will  be  discussed  further  after  some  examples  have  been 
considered. 

Theorem  3 

|x(T,$)|J  <  X'J(0,0)  =  l  .  (18) 

Theorem  4  (Slebert’s  theorem  (4)) 

The  ambiguity  function  is  its  own  two-dimensional  Fourier  transform: 

JJT*  (T.*)|a  **P  t-2wi  drd<t>  =  |x(t./)|3  .  (IS) 

Theorem  2  can  be  obtained  from  Eq.  (19)  with  t  -  f  -  o. 
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Symmetry 

<t>)  -  (  u  (  t )  u*(  f  *  t)  i- iip  ( 2Trii>t)  dt  -  Ju(ij  -  t)  u*(i))  exp  [2  -  r)]  dr) 

=  exp  (  2fff<t>T)  x  *(T'<t’) 

2  •  (20) 

Relationship  to  Matched  Filter 

Consider  a  filter  matched  to  a  real  transmitted  signal  u  ( < );  the  impulse  response 
is  h(f)-  u(T-  t).  Let  the  input  to  this  filter  be  the  time-delayed  and  frequency- shifted 
waveform 

x(f)  =  u(t  i  r  )  exp  f-  2  tt  i  <t>  (  f  ♦  r)]  . 

Then  the  output  is  y  ( t )  -  x ( o  *  />( * )  -  f  x(&)  h(t  -  a.)  da 

y  (  t )  =  fu(iir)  exp  ['2n  i«t>(  a  4  r)]  u  ( T  -  t  +  a)  da 
r  f  u  </{)  u(ft-  t  t  T  -  f  )  exp(-2wj<t>/j)  dp 
V(T)  =  X<t,4>)  . 

Thus,  the  GACF  is  the  output  of  a  filter  matched  to  the  transmitted  signal  in  re¬ 
sponse  to  an  echo  with  constant  time  delay  and  constant  frequency  shift  (narrowband  ap¬ 
proximation  to  Doppler  effect). 

Convolution  Theorems 
Theorem  5 

If  two  functions  are  convolved  in  the  time  (frequency)  domain,  their  generalized 
autocorrelation  functions  are  convolved  in  the  time  (frequency)  coordinate. 

Proof  of  frequency  convolution  case 

Let  u(f)  -  v(f)  »(f).  Then  U(f)  -  V(/)  •»'(/') 

Xu(r,<H)  -  I  V*(  f )  t/(  f  +  <t> )  exp  (  2rrifr)  df 

-  ///  )  *'*(  f-|i)  V(„)  »(/♦  $  -  v)  exp  (  2ri  fr)  dpdvdf 

-  ' J'  -  ^  (>')  If*  (7  )  ftf1/  +  M  -  v  +  <t>)  exp  [2  iti  (>j  +  y  )  t)  d>jdy  dr 

-  ff  *',*(y  )  F(*'  1  V„(7-.y  -  v  t  <t>  )  exp  (  2tm>t)  djjc/i- 

-  )  V(n  4  £)  exp(2wjfir)vw(T,*-  <;  )  dyd£ 

-  I  Xv(T.OV.(r,t-(  )  d< 

=  Vv(T,<t.)  .  V(T,*)  , 
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where  1  indicates  convolution  with  respect  to  <t.  The  proof  of  the  time  convolution  case 
is  similar.  The  operation  of  these  two  theorems  is  noncommutatlve  (31). 


Uniqueness  Theorem 
Theorem  6 


The  function  u(t)  is  uniquely  determined  to  within  a  multiplicative  constant  of  unit 
magnitude  almost  everywhere  (a.e.)  by  its  generalized  autocorrelation  function. 


Proof 

f  Xu(r,Q  )  dr  =  f  I  u  (f  )  u*(  t  -  T )  exp  )  dtdr 

-  i  u(t)  exp(-2iri$t)  (  u*(n)  dqdt  -  U($)U*(0)  . 

Therefore,  if  xu(t.O)  =  xt{r,Q),  V(*)V*(0)  =  0)  and  hence  if  f*(0)  *  0, 

V($)  =  cl/(®),  where  c  =  u*(0)/V*(0).  Since  two  functions  having  the  same  Fourier 
transform  are  equal  a.e.,  we  have  v(t>  =  cu(t),  a.e.,  and 

X'U(T.*)  =  XV(T,«I.)  ^  (  cu  (  ()  c*u*(t  -  t)  exp  (-2"j©t)  dt  =  |  c  |2  Xu(r ,  «> )  . 


so  that  |  c |  s  i. 

Complex  Energy  Density  Function 

Rihaczek  (41)  has  recently  defined  a  complex  energy  density  function  which  may  be 
obtained  as  the  Fourier  tr?msform  of  the  GACF. 

Consider  a  real  signal  represented  by  [u(  t)],  where  u  ( t )  is  an  analytic  signal, 
and  let 


/  |u(t)|J  dt  =  f  |y</)|5  df  =  t  . 

where  £  is  twice  the  total  energy  ci  the  real  signal. 

Then  [  u  ^  t )  |  *  represents  power,  or  "energy  density  waveform,"  and  \U(f)  |s  is  the 
energy  denstty  spectrum.  If  the  autocorrelation  function  of  u(  t)  is  given  by  Eq.  (5),  then 
from  Eq.  (6),  or  by  the  Wiener- Khintchine  theorem,  the  spectral  density  is 

*(')  =  SFt/[*(t)1  =  I^COI2  •  (2D 

Similarly,  the  autocorrelation  function  in  the  frequency  domain,  Eq.  (8),  transforms 
into  the  "waveform  density."  From  Eq.  (9)  we  have 

[«(<*>)}  =  l"(OlJ  • 

By  analogy,  the  two-dimensional  Fourier  transform  of  the  GACF  is 


(22) 
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V  yr/  [V(r.« •)] 


A 


// 


V  exp  I-  2wi(  />-  iff)] 


dr  dt 


-  ///  u  (>/)  u*(  >;  T)  pxp  (-  2»t/4>’7  )  exp  [  2ji/(  /r  -  4>f  )]  dndrdt 
=  u(t)Ju*(t-r)exp(-2rrifr)  dr 
=  u(t)  exp  (-27fi/'0  /  U*(^)  exp  (  2  wi/£) 

=  u  (  t )  U*(  f)  exp  (- 2ni  ft )  .  (23) 

Thus,  just  as  the  autocorrelation  functions  in  time  and  in  frequency  each  are  the 
Fourier  transforms  of  an  energy  density  function,  as  shown  in  Eq.  (21)  and  Eq.  (22),  we 
may  define  the  complex  energy  density  function  in  time  and  frequency  from  Eq.  (23)  as 

-■  u(t)  U*(f)  exp(-2ni  ft)  .  (24) 

Then  the  energy  density  spectrum  is 

/«<*. O  dt  --  |t/(OIJ  ; 

the  energy  density  waveform  is 

/  df  --  |u(t)|J 

and 

)/«(»,/■)  dtdf  -  E  . 
the  total  energy  of  the  analytic  signal. 

The  energy  of  the  analytic  signal  within  a  "cell"  of  area  tb  centered  at  (t0,  f0)  in  the 
time- frequency  plane  is  given  by 

l0*T  2  (q*a  J 

ha  =  /  /  « ( i  •  / )  dfdt 

l0'T  }  /C-B  2 

EXAMPLES  OF  NARROWBAND  AMBIGUITY  FUNCTIONS 
Single  Gaussian  Pulse 

We  consider  first  an  example  discussed  by  Woodward  (1).  Only  slightly  more  general 
than  Woodward's  simplest  pulse  is  the  single  Gaussian  pulse 

u  (  I )  k  i'xp  ( -  .it i)  .  d  0  .  (25) 

where  the  parameter  .-»  determines  the  width  of  the  pulse,  and  k  <  2a  4  so  thci 


I  u2(  t  )  dt  1  . 
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It  Is  shown  in  Appendix  B  that 


.  Xu(T,<t) 


exp  (-nrV 2 


tt2  <t>3/  2a  -  tt 


C 


Thus, 


i\u(T,4.)|2  exp  [-{ arJ  t  u2  <t3/a)]  ,  (! 

and  contours  of  constant  ambiguity  are  given  by  the  ellipse  ar1  +  rr2<t>2/a  =  constant,  as 
shown  in  Fig.  1. 


Fig.  1  -  Curves  of  constant  ambiguity: 


Note  that 


IfluCOl1  -  IVuf-r.O)!5  -  exp(-ar2) 


and 


i"u<«>>!2  =  |X„(0.<t)ia  exp  /a) 

so  that  in  this  case  the  ambiguity  junction  is  factorable:  l\(T.<t>)|2  ;r(t);2  |,  (<t)'2. 

This  result  is  not  general.  When  it  holds  we  have  from  Eqs.  (7),  (10),  and  (17)  that 
TF  l. 


Single  Rectangular  Pulse 

Let  u(t)  -  (2a)’1'3,  ] r |  <  a,  zero  elsewhere.  In  Woodward's  (1)  notation, 
u(t)  -  (2a)’1  3rect(t/2a).  Then 
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Vu(t.4>)  ‘*r ( I )  ’•*<  I  t)]  3t*{(2a)*1  rcct  «2<  -  T)/2(2a-  M)]} 


!(1  *  lT!/2a)|  sine  [( 2.1  -  i  Tj  )  <t]  exp  (-ni<tT),  |  t  |  <2a 

0  .  I  r  |  >  2  a 


I  VU(T.  4>  )  1 2  -  (1  -  |r|/2a)2  sine2  ((2a-  .i  I  <>]  . 


zero  elsewhere,  where  sim-  x  sin  nx. 


Modulated  Uniform  Pulse  Train 


To  facilitate  Woodward's  next  example,  the  Gaussian  pulse  train,  we  consider  a  train 
of  (Dirac)  delta  functions.  Eschewing  questions  of  rigor  and  using  Woodward’s  notation, 
we  let 


w(t)  Y  S(t-nT)  -  repj-S(f) 


W(f)  ^  cunb1/r  1  = 


If  the  pulse  repetition  period  T  is  unity,  the  train  of  delta  functions  Fourier  trans¬ 
forms  into  itself. 

The  GACF  of  a  train  of  delta  functions  is  shown  in  Appendix  C  to  be 


■*>  Z  S(r-nT)  B  -  i). 


This  "bed  of  nails"  (5, 22, 33)  is  shown  in  Fig.  2,  drawn  as  if  the  delta  functions  had  finite 
amplitude. 

Let  a  general  pulse  train  be  represented  by 

v(()  -  rep j-  u  (  t)  t i(  t )  ,  £  5(  t  nT)  -  u  (<)•*(  f)  .  (31) 

By  Theorem  5  we  have 

XV(T.*)  -  X'B(r.4>)*  x-_.(T.<i>) 


:  ^  Z 5  ■  y)  1  a'<t) &  ,r  •  nj~ 

y  Z  *  (•  -■=-)  ■ 


I 
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Fig.  2  -  "Bed  of  Nalls"  generalized 
autocorrelation  function 


If  the  pulse  train  of  Eq.  (31)  la,  in  turn,  modulated  by  an  envelope  function  x(t>,  the 
resulting  waveform  and  Its  spectrum  may  be  represented  by 

y(f)  -  x(t)  v(t)  and  Y  (  f)  =  X  (  f)  .  V  (  f)  . 

Using  Theorem  5  and  Eq.  (32)  it  follows  that 

MT,<t>)  -  Xx(T,®)  *  V,.  (  f ,  ®  ) 

® 

=  /  Xx(r,v)  xu(r-nT,Q-  v)  S  -  ~r  *  ‘'j  dv 

Thus,  the  ambiguity  function  is 


iXy(T,«)| 


it .  .  •,n 


We  now  impose  the  additional  conditions  that  the  GACF  of  a  single  pulse  vanishes 
for  |r|  >  t  and  that  the  GACF  of  the  envelope  function  vanishes  for  |®|  >  l  r.  With 
these  conditions  we  can  write 


V„(T  -tT,®  )  \'*(t  -  nT,<t>) 


|.\u(t  -iT,®)!3 


(33) 


and 


v. 


« 


U  m  ’ 


where  is  the  Kronecker  delta.  Hence, 


NRL  REPORT  7007 


11 


UyCT.V)!1 


(34) 


For  a r -  ctangular  pulse ,  u(f)  -  (2a)‘*  2  rrct  («  2a >,  with  a  rectangular  envelope, 
x ( t )  -  ( 26)' *  1  rect  {t/2b ),  Eq.  (28)  shows  that  the  conditions  Eq.  (32)  are  approximately 
satisfied  if  2a  <  r  «■<  2 b.  We  then  have 


(1-  |t| /2b)a 

72 


£  «t-  |r-nT|/2a)J 

a ,  m 


*  sine 2  (( 2a  -  |t  -  nT|  )  m/T]  si  11c2  [(  2fc  -  |  t|  )  (♦-  m/T)]}  . 


Gaussian  Pulse  Train 

If  u(f)  is  a  single  short  Gaussian  pulse,  we  have  from  Eq.  (26) 
Xu(r-nT,4)xZ(T-mT,*)  =  exp  j-yj^T-n7)a  +  (r-niT)aj  -  4>a  ”  ni$(m-  n)  T 


--  (-1) 


<"->«■  es 


*0 


T)2  +  (T-«r)2 


Z-  *3 


If  n  -  w,  this  equation  reduces  to  I  xu(?  -  nT, «*> )  | 3  in  agreement  with  the  conditions 
in  Eq.  (33).  However,  for  n  »  m,  the  product  does  not  exactly  vanish;  there  is  some 
overlap  of  the  Gaussian  tails.  We  will  have  a  close  approximation  if 


exp 


j-  y  [(T  -nT)2  ♦  (T  -  aiT)2 


«  1  .  Vm  I  fl  .  Vt 


This  requirement  is  satisfied  if  aT2  •  0,  as  must  be  the  case  for  a  short  Gaussian 
pulse  —  moderate  7  requires  large  .3 . 

Similarly,  if  x  (<)  is  a  broad  Gaussian  envelope,  x(»)  k'  «-xp(-/>f2),  we  have  again 
from  Eq.  (26)  with  i>  replacing  a: 

f  3  a 

"  7  <-xp  j-br2  JL  *  ^4>  • 


If 


n ,  this  reduces  to 


If  m  -  n,  we  require  l  bT1  0.  Thus,  for  7  moderate,  we  must  have  h  very  small, 
consistent  with  the  broad  Gaussian  envelope. 

Thus,  if  i>  •  •  7  ■  .1,  the  ambiguity  function  for  a  train  of  narrow  Gaussian  pulses 

with  a  broad  envelope,  obtained  from  Eqs.  (34)  and  (27),  is 
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|xy(T.«t>)|5  £  exp  -a(T-nT)3  ~  ^  (j)  '  bTJ  ‘  V  (^  ’  ^)  '  ^ 

For  each  n  and  m,  the  contours  of  constant  ambiguity  are  ellipses.  Centers  of 
ellipses  are  spaced  at  intervals  of  T  in  r  and  1/T  in  <t>,  and  the  overall  ambiguity  de¬ 
creases  exponentially  in  t  and  $. 

A  sketch  of  this  ambiguity  function  is  given  in  Fig.  19  on  page  122  of  Ref.  1. 

Single  Gaussian  Pulse  with  Linear  Frequency  Modulation 

If  v(f)  =  u(  t)  exp  ( ibt 3)  for  any  u(t),  the  "instantaneous  frequency"  is 

/:  =  i  -  (i)t3)  =  bt/71  . 

*  277  dt  ^  ' 

Thus,  f  i  varies  linearly  with  time.  The  GACF  is 

X„(t.<1>)  =  /v(t)v*(t-r)  exp  (-  277i<J>f)  dt 

-  exp  (-ibr3  )  J u(t)u*(t-T)  exp  [- lvi  (Q-br/v)  t]  d  t 

-  exp  (  ibr3)  X„(T,<t>- br/Tt)  . 

Hence, 

Ixv(t.<D)|s  =  |xu<t,®-6t/7T)|1  .  (?6) 

From  Eq.  (36)  we  see  that 

|Xv(T,br/w)|J  =  |xu(t-.0)|2. 

T.iat  is,  the  ambiguity  along  a  line  $  br/v ,  whose  slope  is  equal  to  the  rate  of  change 
of  the  instantaneous  frequency,  is  equal  to  the  ambiguity  along  the  r  axis  in  the  absence 
of  frequency  modulation. 

For  the  single  Gaussian  pulse  with  linear  FM,  we  have  from  Eqs.  (27)  and  (36) 

ixv(7-,<t>)|3  ‘  exp  j^-sT3  - 

Curves  of  constant  ambiguity  are  again  ellipses, 

fa  +  — )  t3  -  rO  +  II—  4>3  =  constant  ,  (37) 

\  a  /  a  a 

The  eccentricity  of  these  constant  ambiguity  ellipses  depends  on  both  a  and  b,  and 
their  axes  are  rotated  with  respect  to  the  r,<f>  axes  by  an  angle  r>,  where 

tan  2 1)  -  2nb/[v3  -  (a3  *  i>3))  ,  (38) 


as  Is  shown  In  Appendix  D. 
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Note  that  in  this  case 


The  Ideal  Ambiguity  Function  and  its  Approximation 

We  saw  in  Eq.  (27)  and  in  Fig,  1  that  the  parameter  a  suffices  to  reduce  the  ambi¬ 
guity  in  either  the  r  direction  or  the  <t  direction  but  not  in  both.  Equation  (37)  shows  that 
with  two  parameters  it  is  possible  to  reduce  the  ambiguity  both  along  the  r  axis  and  the 
<t  axis  but  not  simultaneously  along  a  line  *  =  r  tan  o.  One  would  surmise  that  with  ad¬ 
ditional  parameters  we  might  do  better.  Ideally,  we  would  like  the  ambiguity  function  to 
be  a  delta  function  at  the  origin,  but  this  we  cannot  achieve.  We  might  hope  to  obtain  as 
an  approximation  a  narrow  spike  in  both  directions  as  shown  in  Fig.  3. 


Fig.  3  -  Approximation  to 
ideal  ambiguity  function 


Equation  (35)  and  Woodward's  Fig.  19  show  that  with  a  train  of  pulses  we  can  reduce 
the  ambiguity  in  all  directions  in  the  vicinity  of  the  origin  to  an  arbitrarily  low  level. 

We  do  this,  however,  at  the  expense  of  having  additional  peaks  appear  elsewhere.  This 
is  a  consequence  of  Theorem  2;  if  we  reduce  the  ambiguity  in  one  place  it  must  pop  up 
elsewhere  so  that  the  total  area  of  the  ambiguity  surface  remains  constant. 

The  closest  to  the  ideal  ambiguity  function  we  can  expect  is  the  so-called  "thumb¬ 
tack"  ambiguity  function.  This  consists  of  a  narrow  spike  surrounded  by  a  uniformly  low 
pedestal,  with  most  of  the  volume  lying  under  the  pedestal.  We  do  not  know  of  any  wave¬ 
form  which  produces  the  "thumbtack"  ambiguity  function.  For  many  applications,  an 
ambiguity  free  region  near  the  origin,  as  obtainable  with  a  train  of  pulses,  is  sufficient. 


Pseudo-Random  Sequences 

We  recall  that  the  ambiguity  function  along  the  t  axis  is  the  square  of  the  autocor¬ 
relation  function  of  the  signal  (see  Eq.  (13)).  Thus,  to  approximate  the  "thumbtack"  am¬ 
biguity  function,  a  necessary  but  not  sufficient  condition  is  that  the  autocorrelation  function 
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of  the  signal  be  small  near' the  origin.  One  practical  signal  known  to  have  such  an  auto¬ 
correlation  function  is  tke  signal  generated  by  maximum  length  pseudo- random  sequences. 
Such  a  sequence,  (n  >,  may  be  specified  by  its  recursion  formula 


cjai-j  •  i=  0.1.2 .  (39) 

7m 

where  cy  =  0  or  1  and  £  indicates  modulo  2  summation  and  by  the  initial  conditions  given 
by  the  values  of 


A  new  sequence  formed  by  the  modulo  2  addition  of  a  maximal- length  pseudo- random 
sequence  and  a  nontrivial  shift  of  itself  will  be  a  shifted  version  of  the  original  sequence. 
That  is,  if 


then 


bj  -  ®  aitT  ,  v«here  T  /  0  (rood  2"  -  1) 

n  n 

=  £  ciai-i  ®  £  ci*uT-i  • 
i-i  i-i 


b ,  -■ 


£ 

i-i 


-j  <*i-j 


®  a. 


► T-J > 


«  t^-J  ■ 

i-i 

Thus,  {bi)  obeys  the  same  recursion  relation  as  <s,).  Since  { «, >  contains  all 
n-tuples  (except  the  all  zero  n- tuple),  these  sequences  are  idential  to  within  a  phase 
shift. 


This  so-called  "shift- and- add"  property  can  be  used  to  obtain  the  narrowband  ambi 
guity  function  of  signals  generated  by  Buch  sequences  for  time  differences  of  r  -  ki0, 
where  *  is  an  integer  and  f 0  is  the  shifting  period  of  the  generator. 

The  additive  group  of  integers  modulo  2  is  isomorphic  to  the  group  consisting  of  -1 
and  1  with  multiplication  as  the  group  operation.  If  we  let  s  ( t )  represent  a  signal  ob¬ 
tained  from  a  pseudo-random  sequence  of  -l's  and  l's,  the  "shift- and- add"  property 
becomes  a  "shift-and- multiply"  property,  as  shown  in  Appendix  E, 

*(t)  s(t-*fp)  -  s(f-  *' t#>  , 

where  *  and  are  Integers.  Thus  from  Eq.  (11) 
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■J 


s  (  t)  s  (t  -  k  t0)  exp  (-  2ni<Pt)  dt 


s  ( t  -  k '  t0  )  exp  (~27,i4>f )  dt 


-  exp  (~2vik  ‘  4>f0  )  S  (C) 


Therefore, 

lx(t<0,4>)|J  -  IS  («t>)|:*  . 

That  is,  for  t  =  ktai  the  narrowband  ambiguity  function  of  maximal- length  pseudo¬ 
random  signals  is  proportional  to  the  power  spectrum  and  is  independent  of  k. 


GENERALIZATIONS  OF  NARROWBAND  AMBIGUITY  FUNCTIONS 
Cross-Ambiguity  Function 

An  extension  of  the  concept  of  ambiguity  functions  to  two  waveforms  was  defined  by 
Stutt  (7). 

Let  Uj(  ()  and  u2<  t)  be  two  complex  waveforms.  We  define  the  generalized  cross 
correlation  function  as 


XjjtT.#)  A  )  u1(t)uJ*(t  -  ' )  exp  (-27ri<J>()  c It 
-  I  t/ *(  /)  Ul(f  t  4> )  exp  (  2n  i  fr )  df  , 

and  the  cross-ambiguity  function  as 

|Xu(T.<t>)|3. 

These  functions  may  be  useful  when  it  is  desired  to  identify  one  of  many  possible  wave¬ 
forms.  Some,  but  not  all,  of  the  properties  of  the  autoambiguity  function  apply  to  the 
cross-ambiguity  function.  For  example,  by  Parseval's  theorem  we  have 

1  !  X„(T,<t>)i3  dll'  -  f  | u, c  *)  I  -  T)|J  dt  . 


(40) 

(41) 

(42) 


so  that 


//  UjjCr.O)!3  ,/<S>di  -  f  iu,(<)i3  /  !Uj(  C  -  T)]3  drdt  -  1  , 
if  the  energy  in  both  waveforms  is  normalized  to  unity.  Furthermore, 

x,j(r.  0)  t  Uj(  r)  UJ(  f  -  T)  dt  . 

\u(0.<t>)  f  I/*(  *  >1)  <11  . 
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but 

|xn(0,o)|J  -  t /  VO^C)  dt\2 

<  f  |ul(0|3  dt  f  |uj(()|3  dt  =  1 

by  the  Cauchy- Schwarz  inequality. 

Therefore,  if  u^t)  *  *u2(t),  where  k  is  a  constant, 

| Xl2 (O.Oil’  <  1  ,  (43) 

in  contrast  with  Eq.  (15).  Also, 

l*jj(T'4’)|a  =  |/  uj(0  uj(  f  -  T)  exp (-2 rriQt)  rft|s 

-  /  !ut(  Ol*  dt  j  |u*(t  -  r)  exp(-2wi<Dt)|2  dt  =  1.  (44) 

It  does  not  follow  from  this  relationship  that  |x'12(t.*)|2  $  |xl2(0,0)]a. 

In  fact,  if  Uj(t)  -  Uj(t-T0)exp(2ni'<ti0t),  t#,  <i>0  *  o,  then  equality  holds  in  Eq.  (44): 
|Xis(To-®o)|J  =  1-  If,  in  addition, 

Uj(t  -  r)  exp  (  2wi<t>f  )  i  *Uj(  t) 

then  from  Eq.  (43), 

|xl2(0,0)|J  <  1  . 

so  that  in  this  case  ix12(0.0) | 2  <  |x,2(V  <t,o>l3  •  In  this  case, 

|Xl2(0,0)|J  *  -r0)exp(2ffi*0f)  uj(t)  dt|J  =  |x„(t# . 4>0 )  | a  . 

By  proof  similar  to  that  for  Theorem  4  it  can  be  shown  that  the  generalization  be¬ 
comes 

//  |X12(t,«)|2  exp  ((-2iri  (/r- ®t>]  drcf®  =  Xu ( t .  /)  X J2 ( * .  t )  .  (45) 

In  place  of  the  symmetry  relation,  Eq.  (20),  we  have 

XjjC-t.-O)  =  /  Uj(  t)  u*(f  *  T)  exp  (2ffi4>t)  dt 

=  exp  (-2ttj<J>t)  /  u,(  t  -  r)  u*(  f )  exp  ( 2v i « t )  dt 
-  exp  (-2t ritr)  X*,(T,4))  . 


Hence 


I  Xl2  (-T.  -® )  |J  =  |x21(r,<f)|J  . 


(46) 
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Convolution  Theorems 
Theorem  7 

Let  Uj(()  =  u2(  t)  =  v2(t)  *  w2<<>.  Then 

Xu1U3(T’4,)  '  **’l''2(T,4>)  . 

7  *  2 

*  x,iVj(t,4>;  .  (47) 

The  proof  is  similar  to  that  of  Th  orem  5. 

Similarly,  if  u; ( t)  =  vt(t)  ( f >  and  u2(t)  v2(t>  »2(f),  then 

X„iUj(T.-J-)  -  XV.v/T.*).  X.^T,*) 

xVj  t, <t> )  *  xWl Vj (t, <t> )  .  (48) 

Invariance  Relations  for  the  Real  and  Imaginary  Parts  of  Ambiguity 
Functions  of  Analytic  Waveforms  (Ref.  7) 

Let  x , j ( t , <t> )  ■  i ij ( t.4> )  +  ifl}{r,< t>)  bethegeneralizedcross-correlationfunctionof 
analytic  waveforms  ut(  o  and  u}(0-  Then 

//  ^JjC7.<l))  drdi>  -  17  f  j2(r,<P)  drd$  -  1'2  ,  (49) 

if  u,(o  and  u2(t)  are  both  normalized. 


Proof 


Since 


II  |v12(t.<Ji)|s  drdQ  -  1 

=  r/  [£J2(t.$)  +  f  ijC1'.®)]  drd*  . 


C,2  =  (l'2)(Xia  +X'*2)  • 

//  CjjCT'f)  drd^  -  (14)//  (vj2  +  2  |Vj2|S  +  X*22)  drd<t> 
-  1/2  +  //  Cx?2  t  X?*)  dTd*  . 


But 
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//  xtj(T,<b)-dTd4>  =  ////  U2(f)V*(f*^)V2(u)U*(v  4$) 

*  exp  [-2 rri  (  f  4  v  )  r]  dfdvdrd^i 


mce  u  (t)  is  analytic,  so  that 


Similarly, 


//  l/}(f)U'(f+  (t)£/2 (-/•)£/*(-/ 4  4>)  dfd<t>  =  0, 


t/  (/)  =  0  ,  /  <  0 


£/.<-/)  =  0  ,  /  >  0 


;/  X.j,(t,<P)  drdQ  -  0  , 


hence, 


J7  <t.4>)  d-rdO  -  i/2  i 


//  €j2(T.$)  =  1/2  ; 


Notice  that  the  autoambiguity  function  is  a  special  case  of  the  cross-ambiguity  func¬ 
tion,  so  that  in  general  the  real  and  imaginary  parts  of  the  ambiguity  function  of  analytic 
signals  contribute  equally  to  the  invariant  volume  under  the  ambiguity  surface. 


The  Most  General  Ambiguity  Function 

Before  considering  other  ambiguity  functions ,  it  is  of  interest  to  reformulate  the 
problem  in  complete  generality.  From  this  generalization  we  will  rederive  the  previous 
results  and  then  define  an  angular  ambiguity  function.  The  generalization  is  used  in  the 
companion  report  to  obtain  a  wideband  ambiguity  function. 

Let  «,(()  and  sa(f)  be  functions,  square  Integra!  s  on  (-®.=°),  representing  signals 
which  we  wish  to  resolve.  Obviously,  if  sx  and  s2  ar<.  to  be  resolved  at  all,  they  must 
differ  in  some  respect.  As  In  Eq.  (1),  we  use  a  mean- square  criterion  to  maximize  their 
difference.  That  is,  we  wish  to  maximize 

d3  =  /  |s,(t)  -  V*)IJ  dt 

=  /  |s,(f  )|3  dt  4  /  |Sj(f)|J  dt  -  2H  [/ dt]  .  (50) 


Equivalently  (as  in  Eq.  2),  we  can  achieve  maximum  resolution  by  minimizing 

iX"j  °r  K-.r- where 

*-x.3  4  f  *}(t)s*^t)dt  . 
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|  v »,  *  2 1 2  may  be  termed  the  "most  general  ambiguity  function."  It  is  too  general  to  be 
useful  in  itself;  its  only  value  lies  in  the  fact  that  from  it  we  can  obtain  less  general  but 
more  useful  ambiguity  functions.  The  argument  and  functional  form  of  Xj,  i2  depend, 
of  course,  on  the  characteristics  of  st  and  x2  that  are,  in  fact,  distinguishable. 

1.  If  ^(t)  -  u(t>  (which  may  be  the  low-frequency  part  of  a  high-frequency  wave¬ 
form  •/>(* )  =  u(  t)  e*p(iwt)  but  need  not  be  so  restricted)  and  if  s3(t)  -  u(t  -  r),  a  time- 
shifted  version  of  the  same  waveform,  we  have 

X.,.,  =  Xu ( r , 0 )  =  lu(t)u\t-r)dt  =  *(r)  . 

the  complex  autocorrelation  function,  Eq.  (5).  In  this  case  |x«,«2|*  =  |*(r)|a,  the  range 
ambiguity  function. 

2.  If  s,(t )  -  u(t )  and  if  s,(  t )  -  u(t )  exp  (2-n  i*4\,  a  frequency-shifted  version  of  the 
same  waveform,  we  have 

X»,.J  =  *u(0.*)  =  /  |u(4)|2  exp  (-  2771*4)  dt  -  *(*)  , 

the  frequency  autocorrelation  function,  Eq.  (8). 

3.  If  Si(t)  -  u(t)  and  if  s2 (4)  =  u(4  -r)  exp  [27ri*(4  -  T)l,  a  time-  and  frequency- 
shifted  version  of  the  same  waveform,  we  have  (neglecting  exp  (  27ri<j>T)) 

X.l<2  =  X^i  .*)  =  J  u(  4)  u*(  4  -  r)  exp  (-  2771*4)  dt  , 


as  in  Eq.  (11). 

4.  If  Xjt  4 )  i  u,(  4)  and  if  ®2(  4)  =  u5(  4  -  t)  exp  [2t7i  *((  -  t>]  ,  a  time-  and  frequency- 
shifted  version  of  a  different  waveform,  we  have 

=  XUlUj(T,*)  -  /  U,(4  )  u*(4  -  r)  exp  (-2771*4  )  <4  4  . 


as  in  Eq.  (40). 


Angular  Ambiguity  Function 

The  distinguishing  features  between  s2  and  s2  need  not  be  temporal.  By  considera¬ 
tion  of  spatial  differences,  Urkowltz  et  al.  (13)  define  an  angular  ambiguity  function. 

We  consider  two  plane  waves  from  distant  sources  incident  on  a  linear  aperture. 

Let  x  be  the  distance  along  the  aperture  from  a  reference  point  (Fig.  4)  and  let  «,  and 
o2  be  the  angles  between  the  directions  o£  arrival  of  the  wavefronts  and  the  normal  to 
the  aperture.  Let  * (()  represent  a  signai  transmitted  with  a  propagation  velocity  c  ,  as¬ 
sumed  constant.  If  the  incide  t  waves  are  echoes  from  stationary  targets  at  the  same 
range, then 
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/ 
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Fig.  4  -  Wavefront  geometry 


Vf>  r  sl(t.x,Oi)  =  s(t  - 
s,(<)  -  Si(t,x,e2)  -  i  - 


x  sin  6. 


x  sin  6. 


and  from  Eq.  (51) 


/  x  sin  0.\  /  x  sin  tf  \ 

X.j.  .0,  )  --  /  *  [t  ~  T - 3 s*  [t  -  r - — -j  dt 

r  *  r  /Sin  ®1  '  Sin  tfa\ 

=  /  set)  S*  \j  ♦  X  ^ - ^jjd 


c/ 1 


If  v  =  ( sin  6j  -  sin  bj  )/e,  then  x,  ,,,(x.  v)  =  f  s(  t )  s*(  t  +  xv)  dt  -  fl^xv),  where  /?,( •  )  is 
the  complex  au'  -correlation  function  of  the  transmitted  signal. 

Urkowitz  et  al.  define  the  angular  ambiguity  function  /(v)  as 


y  (v)  ,  r  |I(x)ia  Ss  (  xv  )  dx  . 


(52) 


where  l(x)  is  the  "illumination  function”  of  the  aperture.  A  detailed  discussion  of  the 
illumination  function  and  its  spatial  Fourier  transform,  the  "space  pattern"  of  the  aper¬ 
ture,  is  beyond  the  scope  of  the  present  treatment. 

If  resolution  in  both  azimuth  and  elevation  is  considered,  a  two-dimensional  illumi¬ 
nation  function  f(x,y)  is  required,  and  a  two-dimensional  angular  ambiguity  function  is 
defined.  If  i(»  v)  l,(x)  I2(y),  the  two-dimensional  angular  ambiguity  function  be¬ 
comes  the  product  of  the  individual  ambiguity  functions. 

If  the  targets  are  at  different  ranges  and  are  moving  with  different  radial  components 
of  velocity,  a  four-dimensional  ambiguity  function  in  azimuth,  elevation,  range,  and  "range 
rate"  (Doppler  shift)  is  required.  For  narrowband  signals  this  four- dimensional  ambi¬ 
guity  function  is  separable  into  the  product  of  the  range-Doppler  ambiguity  function  and 
the  azimuth-elevation  ambiguity  function.  The  same  is  true,  in  general,  of  a  six- 
dimensional  formulation  which  includes  angular  velocity  as  well. 
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If  the  illumination  function  is  complex,  Eq.  (52)  cannot  account  for  the  effect  of  the 
aperture  phase  function  on  angular  resolution.  To  avoid  this  difficulty  and  to  avoid  de¬ 
pendence  of  resolution  on  the  orientation  of  a  receiving  array,  Procopio  et  al.  (15)  replace 
Eq.  (50)  with  the  integrated  squared  difference  criterion 

£2  A  jfl  |  s,  (  t  .0,  .4>j )  -  .  flj  .  ♦j)]2  sin  e  dQdtodt  , 

where 

.«>,)  =  //  I  (x,y)  sl(t,x,y,ei  ,<t>j)  dxdy 

and 

t  ,03  =  J7  l  (x,  y)  s3(  t ,  X,y,e2  ,  *2)  dxdy 

are  the  signals  received  by  the  entire  array,  and  the  angular  orientation  of  the  array  is 
given  by  the  coordinates  o,  <f>.  An  angular  ambiguity  function  may  now  be  defined  as 
iv,u./i>r,  where 

*(<M)  '  ///  S,(f  .»l  <i>i)  <I>3)  sin  0  dddt>dt 


where  a  e1  -  0i  and  ft  -  <S>2  -  ®j. 

Additional  flexibiiity  (and  complexity)  is  obtained  by  allowing  the  illuminating  func¬ 
tion  l(x,y,  f)  to  be  time  varying. 
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Appendix  A 

PROOFS  OF  THEOREMS  1  THROUGH  4 


Theorem  1 

X (7,®)  =  /  u  (  f  )  u*(  t  -  r)  exp  (-2wi$ f  )  df 

=  ///  U(v)U*(f)  exp  {2m  [i  f  -  f  (f  -  t)  -  4>  r  ]  >  drdfdt 
-  If  U(v)  U*(  f)  S[r  -  (  /  v  ®)]  exp  (  2  vifr)  dvdl 
=  JU(f  +  ®)  l/*(  /)  exp (  2m i/V)  df  , 
using  the  well-known  properties  of  the  Dirac  delta  function  £  ( ■). 

Theorem  2 

X(T.*)  -  -  T)J  . 

Hence, 

;  u(t)  U*(f  T)  , 

and  by  Parseval's  theorem, 

/|y(T,®)|3d®  --  /  |u(<)|2  |u(l  -  T)|J  dt  . 

Thu B 

n  |Y  (T,®)  j 3  d®dr  -  f  |  u  (  t  )  |  3  /  |U(T)|3  drdt  -  1  . 

Theorem  3 

| X  c T . 4> ) [ 3  =  1 1  u  ( t )  u*(  t  -  T)  exp  ( -  277 J <t>f  )  dt  ] 1 

£  /|u(f)|3df/|u(f-T)|3  dt  -  1  . 

by  the  Cauchy- Schwarz  inequality. 

Theorem  4 

;ftr*T  t  [l  V(t,®  )|3J  I  J  !y(t,<I>)!3  exp  [-  2m  (  /r  -  ®<  )]  dfd<P 
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^ii ^rt  [!^(t'<|,)|3J  =  ////  "  (1)  u*(1  "  T)  «p(-2wi*n)  u*<0  “  (C ~T) 

x  exp  (2wi*£)  exp  [-  2i7i  (  /V  -  )}  dvd^d-rdS 

-  ff  u(  1)  )  u*(i)  -  T)  u*(n  -  t)  u  (I)  -  t  -  T)  exp(-2 wi/r)  drjd-r 

=  //  u  ( 17 )  u*  ( >)  -  t)u*(£)u(f-  O  exP  [_2wi  ^(>7  -  £)J  dr,d( 
=  |X(t,  0|2  • 


Appendix  B 

GACF  OF  A  SINGLE  GAUSSIAN  PULSE 

u  ( (  )  =  (  2a/v) 1/4  exp  (-at1)  ,  a  >  0 

X  ;  ( la/v)  1  1  /  exp  [-at1  -  a  (t  -  r)2  -  2wi®t]  d t 

-  (2«/e)>'iexp(-aTV2)St#{exp[-2fl(f-r/2),]>  . 

From  Woodward  (1),  page  28,  "pair  3,"  we  have  If,/  [exp  (--nf3)]  =  exp  (-vf1). 
be  shown  as  follows: 

OB 

•f,  /[exp  (-vt 3)]  =  j  exp  (-*( 1  -  2riift )  dt 


-  exp  (-77 


exp  [--nr  (  t  +  i  [  )3]  dt 


at*  i  f 

-  exp(-n/2)  f  cxp(-7rfa)  dt 


-  exp  ( -rr/2 )  I  exp  (~rrt 2)  dt  -  exp  (-nf 2) 

■  Ol 

Then 

if,  /  [exp  (,-2af  3)]  -  (n/2a)'  3exp(-n3/‘3/2a) 

by  Woodward's  "Rule  8,"  and 

If,/  {exp  [-2a(  t  -  t/2)3]  )  -  (7r/2a)'  '  3  exp  [--n1  f1/2a  -  2f‘  /V/2) 


by  "Rule  6."  Applying  this  result  to  Eq.  (Bl)  we  have 

\'(T,4»)  "  exp  ( -ar  3/ 2  -  ti  34>j/2«  -  ttj4>t  )  . 


(Bl) 
Tills  can 
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Appendix  C 

GaCF  OF  A  TRAIN  OF  DELTA  FUNCTIONS 


"(t)  =  £  S(t  -nT) 


X»(T-<t>)  -  f &(t  ~nT)  8  (t  -  t  -  rnT)  exp  (_-2,r’iQt)  dt 


=  2  «[(«-•)  r-T]  exp  (-2»ii.*T)  -  ^  5  (t -nT)  exp  (~2'rrim0T) 


Now  consider  the  formal  Fourier  series  expansion 


y  S(<t>-n/r)  =  2^  Cn  exp  (- 2nJ^<ll7,)  . 


f 

7!  |  5(0- n/T)  exp  (2rrimt>T)  d<P  -  1 


"  -i/ar 


-14 


exp  [-  2tti  (n  -  m  )  OT]  c/t> 


=  £  Cn«n./^  -  Cm/T  , 


so  that  c.  -  T,  Vm.  Hence, 


y  exp  (-2nimt>T)  =  T"1  ^  StO-m/T) 


,(T,0)  ^  7"1  2]]  B  (T  -  nT)  $  (9  -  m/T) 


Appendix  D 


ROTATION  OF  AXES  OF  THE  AMBIGUITY  DIAGRAM  OF  A  SINGLE 
GAUSSIAN  PULSE  WITH  LINEAR  FREQUENCY  MODULATION 


Equation  (37)  will  be  in  the  standard  form  for  an  ellipse  if  we  define  new  axes,  r ' , 
®'  rotated  by  an  angle  e  with  respect  to  the  r,  6  axes,  where 

t  -  t'  cos  9  -  <t>'  sin  6 


and 


-  t'  sin  0  +  <P‘  cos  $  . 


If  this  is  substituted  in  Eq.  (37),  rewritten  in  the  form  At 3  +  b t4>  +  c<t>3  =  K,  and  if 
the  coefficient  of  r'<t‘  is  required  to  vanish,  we  obtain  tan  20  -  B/(A-C)  from  which 
Eq.  (38)  follows. 
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Appendix  E 

MULTIPLICATIVE  PROPERTY  OF  PSEUDO-RANDOM  SIGNALS 


Theorem 

Let  s(o  represent  a  signal  obtained  from  a  maximal- length  pseudo- random  sequence 
of  -Is  and  Is,  of  infinite  duration.  Let  k  be  an  integer,  k  *  o  (mod  p),  where  p  is  the 
period  of  the  pseudo-random  sequence,  and  let  t0  be  the  shifting  period  of  the  generator 
producing  the  signal.  Then  s(t)  s(t  -  ktQ)  =  s(t-k't0),  where  k '  is  also  an  integer. 


»(0  -  E  E-i  o  v«) 

*  i 

Now 


VS2) 


(This  notation  is  similar  to  that  used  on  p.  IS  of  J.  I..  Lawson  and  tl,  E.  Uhlenbeck,  "Threshold  Sig¬ 
nals,”  New  York:  McGrow  Hill, i960.  I  am  grateful  to  Dr.  H.  L.  Saxton  for  calling  niy  attention  to  It. 
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--  aiai-(k-j,i)  =  «,•-*>  •  (E3) 

where  k'  ia  an  Integer  by  the  "shift-and-  multiply"  property  of  pseudo- random  sequences 
of  Is  and  -  ls.t 

Using  Eq.  (El),  we  have 


MO  MO 


1 , 

(( i  -  1 )  fo 

<  t  <  itQ]  n  (< j  - 

O  *0 

0  . 

otherwi se 

1 , 

f0[max  ( i , 

j)  -  1)  <t  <  t0  min 

UJ) 

0  , 

otherwise 

If  i  <  >,  this  product  vanishes  unless  (j  -  l)  f0  <  t  <  it0;  this  means  j  -  1  <  i  <  j, 
which  can  only  fcs  satisfied  for  i  -  j.  Similarly,  if  iz  j,  MO  a,(  t )  =  o  unless 
( i  -  l)  t0  <  t  <  j  f„,  for  which  we  get  i  -  l  <  j  <  i,  and  again  j  -  i. 

Thus,  we  can  write  M  0  MO  =  MOO  >  and  from  Eqs.  (E2)  and  (E3), 


»(O*(t-*t0)  =  E  E  M*«MO«j; 

i  > 


=  E  O-lt'MO 

i 


Q.E.D. 


tSee.  for  example,  C.  McCoy,  Jr.,  '* R>wer  Spectrum  Estimates  of  Sampled  I’seudo-Randoni  Sequences 
NRL  Reixirt  6673,  p.  6H,  Dec.  39,  1967;  Thesis,  The  George  Washington  University,  Feb.  22,  196H. 
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